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ABSTRACT 


C Context. The ESA PLAnetary Transits and Oscillations of stars (PLATO) mission will search for terrestrial planets in the habitable zone of 
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solar-type stars. Because of telemetry limitations, PLATO targets need to be pre-selected. 
Aims. In this paper, we present an all sky catalogue that will be fundamental to selecting the best PLATO fields and the most promising target stars, 


deriving their basic parameters, analysing the instrumental performances, and then planing and optimising follow-up observations. This catalogue 
also represents a valuable resource for the general definition of stellar samples optimised for the search of transiting planets. 
Methods. We used Gaia Data Release 2 (DR2) astrometry and photometry and 3D maps of the local interstellar medium to isolate FGK (V<13) 


and M (V<16) dwarfs and subgiant stars. 


Results. We present the first public release of the all-sky PLATO Input Catalogue (asPIC1.1) containing a total of 2 675 539 stars including 
2 378 177 FGK dwarfs and subgiants and 297 362 M dwarfs. The median distance in our sample is 428 pc for FGK stars and 146 pc for M dwarfs, 
respectively. We derived the reddening of our targets and developed an algorithm to estimate stellar fundamental parameters (Ter, radius, mass) 


from astrometric and photometric measurements. 


Conclusions. We show that the overall (internal+external) uncertainties on the stellar parameter determined in the present study are ~230 K (4%) 
-CÀ for the effective temperatures, ~0.1 Ro (9%) for the stellar radii, and ~0.1 Mo (11%) for the stellar mass. We release a special target list containing 
CO all known planet hosts cross-matched with our catalogue. 
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1. Introduction 


The PLATO mission 2014] [2016) is the third 


medium-class mission in the ESA Cosmic Vision programme. 
Its ambitious goals are the detection and characterisation of ter- 
restrial planets around solar-type stars as well as the study of the 
properties of host stars. The focus of PLATO is on long-orbital- 
period terrestrial planets that are also in the habitable zone of 


oC solar-type stars. PLATO is designed for the discovery of Earth- 
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analogue planets under potentially favourable conditions for the 
development of life. 

PLATO will achieve its goals by detecting the low-amplitude 
dips in stellar brightness produced by planets transiting in front 
of the disk of their parent stars. The mission is designed and op- 
timised to continuously monitor two Long-duration Observation 
Phase (LOP) fields, each one for up to three years, and has the ca- 
pability to monitor some additional fields during shorter times, 
up to three months each (‘Step and stare’ Observation Phase, 
SOP). The precise observing strategy will be decided at a later 


* The catalogue described in this article is only available in elec- 
tronic form at MAST as a High Level Science Product via 
and 
in the SSDC tools page (https:// 
and at the CDS via anonymous 
ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 


u-strasbg. fr/cgi-bin/viz-bin/cat/J/A+A/vol/page 


-mail: Marco.montalto @unipd.1 


stage (two years before launch at the latest), taking advantage of 
the advances in the field at that time. 


PLATO will employ an array of 26 cameras (hereafter ‘cam- 
eras’ signifies the telescope optics and the focal plane assembly, 
including all the ancillary devices like baffling, electronics, etc.) 
made up with refractive 20 cm-class optical elements with a very 
wide field (about 40 degrees in diameter) and a 12cm diameter 
equivalent aperture. 


Each camera field of view is almost covered by an array of 
detectors totaling 1037 square degrees, and 24 cameras are ar- 
ranged in four groups of six cameras each, co-aligned to four 
slightly overlapping fields of view. The S/N of each star will 
therefore be dependent on its precise location within the over- 
all PLATO field of view. The combined field of view of the 
system of cameras amounts to a total of about 2100 square de- 
grees [Magrin et al.2018). The remain- 
ing two cameras are optimised for fast monitoring of very bright 
stars, their colour measurements (e.g. |Grenfell et al.12020), and 
fine guidance and navigation. Their focal plane 1s equipped with 
frame transfer detectors that allow for a coverage somewhat 
larger than 600 degrees. The geometry of the covered patches 
in the sky is kept symmetric, under nominal conditions, for 90 
degree rotations around the line of sight to keep the same cov- 
erage of the field of view throughout a pointing and to allow 
continuous in-flight calibration among different cameras. Com- 
pliance with the performance specifications is computed using 
models for the degradation of the performances over the mission 
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lifetime, mainly due to a moderate degradation of the through- 
put, due to the choice of rad-hardened glasses for the optical 
elements most exposed to cosmic radiation 
Corso et al.|2018), and to the worst scenario in which the detec- 
tors or electrical systems of two of the cameras could fail by the 
End Of Life (EOF). 

Because of limits on the telemetry rate from the satellite to 
Earth, it is not possible to download full CCD images at a suffi- 
cient rate to detect planets; therefore, it is necessary to pre-select 
PLATO targets. For a subsample of targets, customised ‘postage 
stamps’ (imagettes) are downloaded to ground; for the remain- 
ing selected stars, centring and photometry will be performed 
on board and then downloaded to ground. Other, previous space 
missions followed the same approach. The CoRoT Exo-Dat cat- 


alogue (Deleuil et al.|2009), the Kepler Input catalogue (KIC; 
, the K2 Ecliptic Plane Input Catalogue (EPIC; 


, and the TESS Input Catalogue (TIC; Stassun| 
2019) were produced to select the best targets for 
exoplanet searches. 

In this paper, we present the all-sky PLATO Input Catalogue 
(asPIC), whose main purpose is to serve as the primary database 
for the selection of PLATO targets. Establishing the properties 
of PLATO host stars is crucial to determining the properties of 
transiting planets; for example, the uncertainty on their radii is 
directly related to the uncertainty on the host-star radius. The 
classical properties of planet hosts are also essential to modelling 
the stars using asteroseismology. Also, the habitability of a given 
planet depends on the properties of the host star, among other 
factors. Moreover, knowing the values of stellar parameters per- 
mits the statistical study of the frequency of planets as a function 
of stellar type, metallicity, environment, and so on (e.g. [Howard] 
fetal DOT) [Bryson etal. 2020). 

The construction of such a catalogue not only requires the iden- 
tification, across the entire sky, of stars with the required spec- 
tral type and luminosity class, but also requires that noise con- 
straints be accounted for to facilitate the detection of transiting 
planets around the selected sources. Hence limits on the magni- 
tude range of the stellar samples are set according to mission re- 
quirements and the degree of contamination of each target from 
neighbouring sources must be evaluated. The PLATO pixel scale 
will be 15 arcsec pix! on axis, which lies between that of Ke- 
pler (4 arcsec pix ^! ; |Borucki et al.|2010) and that of TESS (21 
arcsec pix^!; . 

The optical quality of the PSF, although slightly variable over 
the field of view and depending on the thermal status of the cam- 
eras, is expected to encompass an area of a few pixels (see also 
Unibriaco et aL DO18). Consequently 
1t is expected that several sources, especially in crowded fields, 
will be blended into a single brightness measurement area and 
they may become the cause of false positive signals (e.g. 


2015 2013). More specifically, the as- 
PIC will be used to: 


1. select the optimal PLATO observing fields; 

2. select all FGKM dwarf and subgiant stars satisfying the mag- 
nitude and signal-to-noise constraints defined in the PLATO 
stellar sample requirements; 

3. estimate basic stellar parameters for all targets, such as tem- 
peratures, radii, and masses; 

4. identify known variable stars, binaries, members of multiple 
systems, active stars with bitmasks, and summarize the in- 
formation available from existing catalogues, that is, mainly 
Gaia and TESS releases; 

5. supply a list of all target contaminants up to a specific angu- 
lar distance from the target and up to a limiting magnitude, 
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for the targets in the fields to be observed, in order to under- 
stand their impact on photometric performances and because 
of the challenge that background contaminants pose to plan- 
etary candidate validation; 


6. guide the organisation and optimisation of ground-based 
follow-up strategies. 


Previous attempts to produce this catalogue showed that this is 
not a trivial task. In particular, it is hard to robustly discriminate 
between dwarfs and evolved stars without accurate knowledge 
of the absolute luminosities of sources. Reduced proper motions 


were employed in the past to this purpose (Nascimbeni et al. 
2016), but this technique is now superseded by the availability 


of trigonometric parallaxes from Gaia DR2 
essentially for all targets of PLATO interest, 
permitting accurate determination of distances and therefore of 
absolute luminosities. 


After ESA adoption (in June 2017), PLATO recently entered 
the full industrial development phases C-D (March 2020) to be 
ready for launch in 2026. At this development stage, the selec- 
tion of the PLATO LOP fields is not yet finalised (it will be two 
years before launch). The present work is focused on the asPIC. 
This catalogue is made publicly available to the community and 
it will be updated whenever relevant changes occur, including 
new Gaia releases. 


After briefly describing the mission requirements defining 
PLATO stellar samples in Sect. 2] we describe the target selec- 
tion criteria in Sect. [3] Then, in Sect [4] we describe how we ac- 
count for reddening and absorption in our analysis. Our stellar 
parameter pipeline is presented in Sect.[] In Sect|6] we compare 
our results with those published in various literature sources. In 
Sect.[7] we describe a special list of known planet host stars in- 
cluded in the catalogue release. Section|8|provides a brief sum- 
mary of the present work. 


2. PLATO stellar sample requirements 


The four main PLATO stellar samples (named P1, P2, P4, and 
PSP]win be observed in the wavelength range between 500 nm 
and 1000 nm, and an additional colour sample will be observed 
in two broad blue and red spectral bands spanning the 500-675 
and 675-1000 nm wavelength ranges, respectively 
2018). The science requirements for the PLATO stellar sam- 
ples are listed below and summarised in Table |1| according to 
the ESA Science Requirements Document (SciRD, PTO-EST- 
SCI-RS-0150_SciRD_7_0). The definition of the PLATO sam- 
ples requires knowledge of the visual apparent magnitude V be- 
cause the SciRD adopts this magnitude as reference to identify 
sufficiently bright targets for ground-based spectroscopic follow- 
up. The V magnitude we used to select the targets comes from 
our own calibrated transformation of Gaia DR2 magnitudes and 
colours to the Johnson V band (Eq. p. This choice ensures the 
maximum homogeneity of the V band photometry. The calibra- 
tion procedure is described in Appendix[A] The colour transfor- 
mation relation is 


! Although Gaia EDR3 has now been released (Gaia Collaboration 


this work is based on Gaia DR2. Future updates of the 
asPIC will be based on Gaia DR3. 

? For historical reasons, the sample P3 has been eliminated, but the 
numbering of the PLATO samples was left unchanged. 
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Table 1. Summary of science requirements for the PLATO stellar samples. 


P1 P2 P4 P5 Colour sample 
Stars 215 000 (goal 20 000) 21000 25000 2245 000 300 
Spectral Type Dwarf and Dwarf and M Dwarfs Dwarf and Anywhere in 
subgiants F5-K7 subgiants F5-K7 subgiants F5-late K the HR diagram 
Limit V 11 8.5 16 13 - 
Random noise (ppm in 1 hr) «50 «50 Š 2 = 
Wavelength (nm) 500-1000 500-1000 500-1000 500-1000 Blue (500—675) 
red m5 o0 
spectral bands 


(G — V)o = co + ci(Gnp — Gnp)o + Co(Gpp — Gnp)o* 
+ ¢3(Gpp — GrP) + c4(Gap — Gnp)o* 


+ c5(Gpp — Gnp)o. 


a) 


where G, Ggp, Grp are the Gaia DR2 magnitudes and the null 
subscript indicates colours dereddenend with the procedure de- 
scribed in Sect.[4] The best-fit coefficients of the transformation 
relation are reported in Table 2] 


2.1. Stellar sample 1 (P1) 


— The total number of targets in stellar sample 1 (cumulative 
over all sky fields) shall be at least 15 000 dwarf and subgiant 
stars of spectral types from F5 to K7, with a goal of 20 000. 


— The dynamic range of stellar sample | shall be V < 11. 


— The random noise level for stellar sample 1 shall be below 
50 ppm in 1 h. 


— In stellar sample 1, the proportion of brighter targets (V < 
10.5) shall be maximised. 


— Stellar sample 1 shall be observed during a LOP. 


2.2. Stellar sample 2 (P2) 


— The total number of targets in stellar sample 2 (cumulative 
over all sky fields) shall be at least 1000 dwarf and subgiant 
stars of spectral types from F5 to K7. 


— The dynamic range of stellar sample 2 shall be V < 8.5. 


— The random noise level for stellar sample 2 shall be below 
50 ppm in 1 h. 


— Stellar sample 2 shall be observed during a LOP. 


2.3. Stellar sample 4 (P4) 


— The total number of targets in stellar sample 4 (cumulative 
over all sky fields) shall be at least 5000 cool late-type dwarfs 
(M dwarfs) monitored during a LOP. 


— The dynamic range of stellar sample 4 shall be V < 16. 


2.4. Stellar sample 5 (P5) 


— The total number of targets in stellar sample 5 (cumulative 
over all sky fields) shall be at least 245 000 dwarf and sub- 


giant stars of spectral types from F5 to late K. 
— The dynamic range of stellar sample 5 shall be V < 13. 
— Stellar sample 5 shall be observed during a LOP. 


2.5. Colour sample 


— Part of the payload must provide photometry in at least two 
separate colour broad-bands. 

— PLATO shall have the capability to observe 300 stars (in 
at least two pointings) located anywhere on the HR dia- 
gram, with colour information, imagettes, and a sampling 
time equal to 2.5 s. 


The above stellar samples are defined only across the two 
selected PLATO LOP fields. The final decision on these fields 
will be taken two years before launch. So far, two preliminary 
LOP fields have been identified. The P1 sample is a subsample of 
P5. However, the P1 sample is also S/N-limited and, because of 
the particular configuration of PLATO pointing, it can be defined 
only once the LOP fields are selected. 

In this work, we construct an all-sky catalogue, focusing on 
only two main PLATO stellar samples. One sample is related 
to FGK (limited to F5) dwarf and subgiant stars with apparent 
visual magnitude V < 13 (and formally corresponds to the P5 
sample as defined above), while the other sample is related to 
M-dwarfs with V < 16 (which corresponds to P4). For the sake 
of clarity, hereafter we refer to these two samples as the FGK 
sample and the M sample, respectively. We note that although 
the asPIC is an all-sky catalogue, PLATO will not point towards 
the Galactic plane. 


3. Selection criteria 


The selection of the PLATO targets was performed using the ab- 
solute colour-magnitude diagram (CMD) corrected for redden- 
ing. First, we employed theoretical models based on Galactic 
simulations to identify the region of the diagram where our tar- 
gets are located (Sect.B-1). On the same diagram, we then repre- 
sented stars that, based on available spectroscopic and/or astro- 
metric observations, have the properties required by the SciRD 
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Table 2. Coefficients of the polynomial best-fit model relating the (G — V)o and (Ggp-Ggp)o colours and residuals of the fit (see also Appendix[A]. 


The relation is valid within 0.5«(Ggp-Ggp)o «5. 


(Gpp-Grp)o (GaP-Gnp), (Gmr-GnpPy, (GarP-Gmp); (Gmp-GRP) Crh 
co Cl C2 C3 C4 C5 
(G-V)y  -0.17276 0.47885 -0.71953 0.24374 -0.04458 0.00317 0.02745 


(Sect. [3.2). Finally, by taking into account both the theoretical 
and empirical distributions of stars in the CMD, we proceeded 
by defining an analytical selection region, and estimated its de- 
gree of completeness and contamination (Sect. |3.3). Finally, in 
Sect. 3.5] we show that our selection does not imply any bias in 
the metallicity distribution of the selected stars with respect to a 
sample of local stars (Sect. [3.4) and present the asPIC1.1 cata- 
logue resulting from our selection. 

As summarised in Sect. |2| the PLATO samples considered 
here include an FGK dwarf and subgiant sample (limited to 
F5), and an M dwarf sample. Considering the definitions re- 
ported in Table 5 of (2013), spectral type 
F5 corresponds to Tef = 6510 K, or to an unreddened Gaia 
colour (Ggp-Ggp)o = 0.587, and spectral type MO corresponds 
to Ter = 3870 K or unreddened (Ggp-Gngp)o = 1.84, where Ggp 
and Grp are the blue and red Gaia magnitudes, respectively. In 
addition, dwarfs and subgiants have log g > 3.5. 


3.1. Theoretical models 


To determine and calibrate our selection criteria in an ob- 
servational CMD, we used Galactic simulations from TRILE- 
GALvI.6, TRILEGAL is a population 
synthesis code for simulating the stellar photometry of any 
Galaxy field. It allows the user to simulate the photometry in 
several different broad-band photometric systems. Each simu- 
lation we performed corresponds to an area of ten square de- 
grees. The simulations were centred at (b)=(65°, 30°) and 
(1,b)=(253°, -30°), and on four surrounding locations ten de- 
grees apart from the central fields in Galactic longitude or lat- 
itude. The two central fields correspond to the two provisional 
PLATO long-duration fields in the Northern and Southern Galac- 
tic hemispheres, that is the North PLATO field and the South 
PLATO field (NPF and SPF, respectively). We considered the 
Gaia DR2 passbands of [Maíz Apellániz & Weiler| (2018) and 
determined the apparent visual magnitude V using the relation 
reported in Appendix [A] We corrected the apparent magnitudes 
of each simulated star for extinction using the V band extinction 
tabulated in the simulations and the relations given in Sect. [5.2] 
in order to convert the extinction in the V band to the extinction 
in the Gaia bands. Figure [I| shows the intrinsic colour (redden- 
ing free) and absolute magnitude obtained from TRILEGALv1.6 
simulations. The left panel shows stars with V<13, log g>3.5 and 
3870 K < Teg < 6510 K corresponding to our adopted limits for 
the FGK sample, while the right panel shows stars with V<16, 
log g>3.5 and Tef < 3870 K (M sample). Figures 2]and B|show 
the distributions of effective temperatures (left panel) radii (mid- 
dle) and masses (right) of the simulated stars. Magnitude limited 
samples, like the ones we are considering, are dominated by hot- 
dwarfs and subgiants, because these stars are intrinsically more 
luminous and therefore observable in a large volume. 

We also repeated the simulations for fields located in the 
Galactic plane or at the Galactic poles. The intrinsic colours and 
magnitudes of the selected stars are perfectly compatible with 


? http://stev.oapd.inaf.it/cgi-bin/trilegal 
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the ones shown in Fi gure[1] Larger uncertainties in the reddening 
estimates occur for the Galactic plane and would lead to higher 
levels of contamination, particularly for hot dwarfs and subgiants 
(see also Sect. [3.3]. 


3.2. Empirical control samples 


For FGK dwarfs and subgiants, as an empirical control sample, 
we used stars from the RAVE DRS survey (Kunder et al.|2017). 
RAVE provides two sets of stellar parameters: those that are di- 
rect outputs of the stellar parameter pipeline and the calibrated 
parameters. Because the direct outputs are based on a grid of 
synthetic spectra which may not match real spectra over the en- 
tire domain of the parameter space, a calibration based on refer- 
ence stars is performed in order to minimise systematic offsets. 
For example, surface gravities are calibrated using both the as- 
teroseismic log g values of 72 giants from Valentini et al.|(2016) 
and the Gaia benchmark dwarfs and giants (Heiter et al.[2015). 
We then adopted the calibrated parameters. We considered stars 
with RAVE parameter ALGOCONV=0 and S/N > 50 to ensure 
the good quality of extracted parameters. A total of 122 167 stars 
were used. We analysed the distribution of RAVE stars in the ab- 
solute colour magnitude diagram constructed using 
distances. We corrected the apparent magnitudes for 
extinction using the procedure presented in Sect. To identify 
the selection regions of dwarfs and subgiants in the (Gpp — Gpp)o, 
Mg plane, we considered a grid of 0.05 mag x 0.5 mag bins in 
colour and magnitude across the plane. We considered stars with 
V«13, where the V band magnitude was calculated with the cal- 
ibration relationship reported in Appendix Stars for which 
RAVE parameters did not satisfy the conditions log g> 3.5 and 
3870 K< Teg «6510 K were considered contaminants. We then 
selected those grid rectangles with at least 50 stars and for which 
the number of dwarfs and subgiants was larger than 75% of the 
total number of stars. Figure |4| (left panel) shows the result of 
this procedure. The black rectangles are the regions belonging to 
the dwarf--subgiant class. 

For M-dwarfs, we used The REsearch Consortium On 
Nearby Stars (RECONS?) catalogue (Henry et al. [2018], and the 
catalogue. Among other information, the RE- 
CONS catalogue reports the apparent visual Johnson magnitude 
of the targeted stars. Whenever multiple measurements for the 
same object were present in the catalogue, we averaged their vi- 
sual magnitudes. From the initial list of 348 entries we obtained 
a list of 289 measurements for the individual objects. The cata- 
logue was cross-matched with Gaia DR2 and we then selected 
stars with V < 16 that satisfy Eq. (1) and Eq. (2) of[Arenou et al.] 
(2018) Gaia astrometric and photometric quality criteria and are 
classified as M-dwarfq*| restricting the list to 195 stars. We then 
calculated the reddening as described in Sect.|4] The right panel 
of Figure [4] shows the absolute colour-magnitude diagram cor- 


^ We also compared our reddening E(B-V) with the one calculated in 
RAVE and found perfectly consistent results: AE(B-V)=(-0.01+0.02). 

5 www.recons.org 

$ RECONS does not report the effective temperature, but it reports the 
spectral type of the stars. 
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9 -  V«16 AND log g>3.5 AND T 


str £3870 7 


V<13 AND log g>3.5 AND 3870 < T, «6510 


0.5 1 1.5 2 2 3 4 5 
(Ga — Grp)o (Gap 7 Gap)o 


Fig. 1. Absolute magnitude vs. reddening-free colour obtained from TRILEGALv1.6 simulations. Crosses represent simulated stars from a total 
area of 100 square degrees in the NPF and the SPF, satisfying the conditions described by the labels at the top of each plot. The left and right panels 
also show the analytical selection (red lines) described in Sect. [3.3] for FGK dwarfs and subgiants (FGK sample) and for M dwarfs (M sample), 
respectively. In the diagram on the right, the blue dashed line represents a 10 Myr solar metallicity isochrone from the Padova database 


fet al.[2012). 
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Fig. 2. Distribution of effective temperatures (left panel), stellar radii (middle), and masses (right) for FGK sample stars from TRILEGAL simu- 
lations. 
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Fig. 3. Distribution of effective temperatures (left panel), stellar radii (middle), and masses (right) for M sample stars from TRILEGAL simulations. 
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rected for reddening of the selected stars (open black circles) of 
this sample. 

constructed a sample of nearby late K and 
M dwarfs (all within 40 pc from the Sun) for which optical and 
near-infrared spectra were obtained. We used the effective tem- 
peratures derived by these latter authors to isolate only M dwarfs 
by imposing Teg < 3870 K. Gaia DR2 was cross-matched with 
this catalogue and we restricted the sample to those stars satis- 
fying Eq. (1) and Eq. (2) of (2018). The apparent 
visual magnitude in the Johnson system (necessary for our se- 
lection, and not reported in the catalogue) was calculated as de- 
scribed in Appendix [A] By means of these selection criteria the 
initial sample of 179 stars was reduced to 153 stars represented 
by the magenta asterisks in Figure [f] 


3.3. Analytical selection 


For the selection of M and FGK targets from the Gaia CMD, we 
then proceeded by defining a set of inequations. The red colour 
limit corresponding to (Ggp-Grp)p=1.84 was used to separate 
M-dwarfs from FGK dwarfs and subgiants following [Pecaut &| 
[Mamajek](2013), while the blue colour limit of the FGK samples 


was set at (Ggp-Grp)o=0.56 corresponding to a spectral type be- 


tween F4V-F5V according to Pecaut & Mamajek| (2013). This 
limit is bluer than the (Ggp-Grp)9=0.587 proposed by |Pecaut 


for F5 stars in order to (partially) account 
for reddening errors. By using TRILEGAL simulations, we per- 
turbed the reddening correction according to the expected red- 
dening uncertainties derived from our map (see Sect. [4), and es- 
timated the number of contaminants (stars with spectral type ear- 
lier than F5) and targets in narrow colour intervals of 0.05 mag. 
Because our samples are magnitude limited, hotter stars are more 
numerous than cool ones, because they are intrinsically more lu- 
minous (see e.g. Fig.D] left panel). The major source of contami- 
nation comes from massive dwarfs and subgiants of spectral type 
earlier than F5. Figure [5] (top panel) shows a colour—magnitude 
diagram from Galactic simulations presenting dwarfs and sub- 
giant stars with V<13 and Tef «7220 K (i.e. later than FO). At the 
limit of (Ggp — Grp)o=0.56 (dashed red line) simulations show 
that we expect a loss of ~10 % of good targets, but an inclu- 
sion of ~90 % of contaminating stars. Extending the selection 
criteria to bluer colours would increase the number of contam- 
inating stars (i.e. stars with spectral type earlier than F5), with 
negligible recovery of good targets (bottom panel of Fig. [5). The 
upper luminosity threshold of FGK subgiants was set to include 
stars down to log g=3.5. Such a definition also helps to reduce 
the contamination from evolved, reddened sources. The lower 
luminosity threshold of FGK dwarfs is set at about three magni- 
tudes below the main sequence to avoid contamination from po- 
tentially spurious sources populating regions of the CMD where 
no plausible main sequence dwarf is expected to be found. With 
our selections, distant red giants are excluded by construction 
independently of their reddening and the same holds for white 
dwarfs. 

In general, theoretical predictions and observational data 
agree reasonably well. For M-dwarfs, such as those visible in 
Figs.[1]and[4] simulated stars appear systematically offset toward 
larger absolute magnitudes than observations, and appear to span 
a smaller colour range. Some differences are also evident among 
the observational samples. Several RECONS stars are intrinsi- 
cally brighter than the[Mann et al.|(2015) sample stars. The upper 
absolute luminosity limit of the M sample was set by perform- 
ing a linear regression of a 10 Myr solar metallicity isochrone 


from the Padova database (Bressan et al.[2012), represented by 
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the blue dashed line in Fig. |1| Considering the theoretical and 
observational uncertainties, this threshold is a good compromise 
to include all cool late-type dwarfs at any metallicity, including 
binaries and active stars. At the same time, by adopting this limit 
we are more likely to discard giant stars, which are expected to 
be brigher than this limit. To reduce the contamination from red- 
dened stars lurking inside our selection region we also limited 
the distance of the stars in the M sample to 600 pc. The rationale 
behind this choice is related to our limiting magnitude V=16. 
Assuming the conversion between Johnson V and Gaia magni- 
tudes in Appendix [A] the maximum distance at which we can 
detect an unreddened MO type star with an apparent visual mag- 
nitude V=16 and at the bright limit of our selection (Mg 96.55) 
is 573 pc. Therefore, any other M-type star within our selection 
region cannot be found at a distance of greater than 600 pc. In 
summary, our adopted selection criteria are defined as follows: 


(Ggp — Grp)o = 1.84 
Meo > 2.334 (Ggp — Grp)o + 2.259 


Mees Distance < 600 pc , (2) 
V < 16 
0.56 < (Gap — Gnp)o < 1.84 
FGK sample = Mcp € 4.1 (Gap — Ggp)o + 5.0 (3) 


Moco > 4.1 (Ggp — Ggp)o — 2.2 
V<13 


The quantities that appear in the relations above are the intrinsic 
colour (Grp-Gpp)o in the Gaia bands, the Gaia intrinsic absolute 
magnitude (Mg), the apparent visual magnitude (V), and the 
distance from [Bailer-Jones et al.| (2018). The intrinsic colour is 
determined by deredenning as discussed in Sect.|4| The selection 
regions are illustrated in Fig. [I] red lines) and in Fig. [4] (dashed 
lines). 

By adopting this analytic selection we find that 92% of the stars 
whose RAVE parameters are compatible with the selection crite- 
ria are included. For the M-dwarf samples we find that 97% (190 
out of 195 stars) of the stars, and 99% (152 
out of 153 stars) of the (2015) stars are included in 
the selection. 

An alternative estimate of the completeness and contamina- 
tion can be obtained using theoretical models. We used the same 
simulations described above in this section to determine the true 
positive rate (TPR) and the false positive rate (FPR) implied by 
our selection. We distinguished the cases of M and FGK selected 
stars. A star is considered a true positive (TP) if it satisfies in- 
equation []for the M sample or inequation3|for the FGK sample 
and the corresponding theoretical constraints ( V<16, log g>3.5, 
and Teg < 3870 K for M sample and V<13, log g>3.5, and 3870 
K < Teg < 6510 K for FGK sample); it is a false negative (FN) if 
it does not satisfy inequation|2] (or inequation[) but does satisfy 
the theoretical constraints; it is a false positive (FP) if inequa- 
tion [2|(or inequation[3) is satisfied but the theoretical constraints 
are not; and finally, it is a true negative (TN) if both inequationp] 
(or inequation B) and the theoretical constraints are not satis- 
fied. The TPR and FPR are therefore defined as TPR= Z5; and 


FPR= PN The TPR and the FPR are measures of the perfor- 
mances of the selection algorithm. The TPR of the FGK selec- 
tion is 100% with a FPR of 12%. The TPR of the M selection is 
88% with a FPR of 0%. The FGK selection appears more ‘per- 


missive’ than the M selection, implying higher completeness but 
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Fig. 4. In the left panel, black rectangles represent FGK dwarfs and subgiants from RAVE (see Sect. [3.2] 
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shows M-dwarfs (Henry et al.]2018| open circles) and (2015| magenta asterisks) catalogues. Dashed lines in both panels show the 


analytical selection defined in Sect. for the FGK (left) and M stellar samples (right). 


also higher contamination. Most of the decrement in the TPR of 
the M sample appears to be related to the sharp dependence on 
metallicity of stellar models of late-type stars. 


3.4. Impact of selection on the metallicity distribution 


It is important to understand if the selection criteria we adopted 
introduced any appreciable bias on the metallicity distribution 
of the target stars. To this purpose, we retrieved the Geneva- 
Copenhagen Survey (GCS) catalogue and 
cross-matched it with asPIC by angular separation (accounting 
for proper motions) and brightness (by comparing the V magni- 
tude of GCS with the V magnitude derived from Gaia; see Ap- 
pendix [A). In particular, we considered two samples, one with 
all stars in the GCS having Teg « 6500 K and for which the 
value of the metallicity [Fe/H] was reported (in total 12 652 
stars), and another one with all stars in the GCS for which the 
asPIC selection criteria (Eq. p] and Eq. |3) were satisfied (and 
had a metallicity estimate; in total 10 798 stars). We then com- 
pared the corresponding distributions of metallicities as shown in 
Fig. [6] We performed a Kolmogorov-Smirnov two-sided test to 
evaluate whether or not the two distributions were drawn from 
the same population. We consider a significance level equal to 
a =0.01 and obtain a p-value equal to 0.6625. The null hypoth- 
esis is therefore not rejected and the two distributions are statis- 
tically equivalent. This allows us to conclude that our selection 
criteria are not biasing the distribution of metallicities of the tar- 
get stars (with respect to the distribution of metallicities of the 
GCS catalogue). 


3.5. asPIC1.1 


Figure [7] illustrates the selection in the Mg versus (Ggp-Ggp)o 
diagram while Fig. |8| (left panel) presents the distributions of 
distances for the M and FGK samples. The median values of 


M and FGK distances are 146 pc and 428 pc, respectively. The 
right panel of Fig. [8] shows the distribution of relative distance 
errors for the different stellar samples. The error on the dis- 
tance is calculated from the upper and lower distance limits re- 
ported by and is equal to the semi- 
difference of these values. In terms of percentage errors in dis- 
tance, the median values for the M and FGK stars are 0.6 and 
1.696, respectively. Considering the median distances of the stel- 
lar samples reported above, this means that M and FGK dis- 
tances are set with a precision of 0.9 pc and 6.8 pc, respec- 
tively. Figure p] shows the distribution of the selected FGKM 
dwarf and subgiant stars across the celestial sphere in a Galac- 
tic coordinate reference system and the position of the two 
currently provisionally defined LOP fields 
(2016). These fields are centred at Galactic coordinates 1=65°, 
b=30° (RA=17":40™:19.265°;, DEC=+39°:35’:01.47”) for the 
NPF and 1=253°, b=-30° for the SPF (RA=05":47™:11.700°; 
DEC=-46° :23’:45.37’”). The definitive choice of the LOP fields, 
to be frozen at the latest two years before launch, involves a com- 
plex optimisation task to merge several constraints and priorities 
of both engineering and scientific nature and will be discussed 
in Nascimbeni et al. (in preparation). Nevertheless, the two pro- 
visional fields illustrated in Fig. [9]already satisfy all PLATO re- 
quirements (Sect. 2). 


The final all-sky catalogue is named asPIC1.1 and contains 
2 675 539 stars. There are 2 378 177 FGK dwarfs and sub- 
giants (V <13) and 297 362 M-dwarfs (V <16). In Appendix[B] 
we describe an alternative selection of FGK dwarf and subgiant 
stars based on Johnson B and V photometry, which was used 
to construct the first version of this catalogue (asPIC1.0). In 
Appendix |C| we describe the catalogue content, while in Ap- 
pendix [D] we present the details of the implementation proce- 
dure that has been adopted to construct asPIC1.1. 
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Fig. 5. Top: Colour-magnitude diagram from Galactic simulations pre- 
senting dwarfs and subgiant stars with V<13 and Tez «7220 K (i.e. 
later than FO). The arrow indicates the average direction of the redden- 
ing vector. Bottom: Percentage of targets (continuous line) and contam- 
inants (i.e. stars with spectral type earlier than F5, dotted line) within 
0.05 mag wide colour intervals. The red dashed vertical lines denote 
our adopted blue colour limit of the PLATO FGK sample, the blue dot- 
ted vertical lines indicate the colour where the selection is 100% dom- 
inated by contaminants, while the blue long-dashed vertical lines show 
the colour where the selection is 10096 dominated by the targets. 
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4. Reddening and absorption 


Because of the presence of interstellar matter, stars appear redder 
and fainter than they are. This has an impact on the selection of 
PLATO targets which is based on the use of photometric quan- 
tities like colours and magnitudes, as discussed in Sect. In 
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Fig. 6. Metallicity distribution of stars selected from the GCS catalogue 
(Holmberg et al.|2009). The black histogram shows all stars from the 
catalogue (with Ter «6500 K) while the red-dashed histogram shows 
all stars selected according to asPIC criteria (inequation D]and inequa- 


tion[3). 
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Fig. 7. Combined selection of FGK and M stars in asPIC1.1. 
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particular, as the PLATO samples are magnitude limited, F-type 
stars, being intrinsically more luminous are located at larger dis- 
tances on average than later type stars, and are therefore the most 
affected by reddening. By neglecting the reddening and absorp- 
tion correction we estimate that 100% of the stars at the blue 
limit of our colour selection (Grp-Gpp=0.56, Eq. would be 
contaminant dwarfs and subgiants (i.e. hotter than spectral type 
F5). It is therefore important to account for the impact of inter- 
stellar matter before selecting the target stars. 
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Fig. 8. Left: Distributions of distances for M (red histogram) and FGK (blue histogram) stars as asPIC1.1. Right: Distribution of relative errors on 


distances in asPIC1.1. 
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Fig. 9. Distribution of asPIC1.1 dwarf and subgiant FGKM stars across the celestial sphere (grey points) in a Galactic coordinates reference 
system. The plot also shows the position of the two provisional long-duration fields (NPF and SPF). 


Different approaches can be considered to determine redden- 
ing and extinction. One of the drivers for the choice of the best 
approach to follow for the PLATO input catalogue is to consider 
that the most important targets lie in the solar neighbourhood. 
Considering FGK targets, reddening should be estimated up to a 
distance of about 1 kpc. 

Figure[10|shows the cumulative distributions of reddening (left) 
and HIPPARCOS distances (right) for all stars in the Geneva- 
Copenhagen Survey (GCS]/] The figure shows that 90% of these 
stars are within 150 pc and the reddening is E(b — y) x 0.02 cor- 


7 We used here the GCSII version (Holmberg et al.|2007) which reports 


the reddening. 


responding to E(B — V) < 0.04. In the GCS, reddening is es- 
timated from the calibration by [Olsen] of the Strómgren 
photometry intrinsic colour index (b — y). The stated precision is 
0.009 mag. The same figure also shows that around 2046 of the 
stars have negative reddening, and are too close to allow a reli- 
able reddening determination At the small distances of PLATO 
targets, determining the reddening on a purely empirical photo- 
metric basis is challenging and requires measurements of very 
high precision. Such estimates can be obtained from measure- 
ment of interstellar neutral sodium absorption lines imprinted 
on spectra of background stars (e.g. [Vergely et al./2001). Such 
measurements are best performed on early-type stars and extrap- 
olation techniques can be used to reconstruct tomographic maps 
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Fig. 10. Left: cumulative distribution of reddening values in GCS. 
Right: cumulative distribution of HIPPARCOS distances in GCS. Hori- 
zontal lines denotes the 90" percentile of the cumulative distributions. 


of the local interstellar medium in all directions. Such spec- 
troscopic measurements have a superior precision than purely 
photometric estimates, but the empirical database is still limited 
(Welsh et al.[2010), and the resulting reddening maps have a lim- 
ited spatial extent (« 300 pc from the Sun) and would be largely 
incomplete for our samples. 


In general, spectroscopic measurements can be used to deter- 
mine the effective temperature and therefore the intrinsic colours 
of stars, which can then be compared with the observed colours 
in several optical and near-infrared bands to infer colour ex- 
cesses (Chen et al./2019). These methods suffer from the com- 
pleteness problem mentioned above (although the availability of 
large spectroscopic surveys mitigates this problem), from the in- 
homogeneity of the spectroscopic samples adopted and from the 
different accuracies and precisions of the effective temperature 
determinations. 


A reasonable compromise among the different techniques can 
be obtained by coupling empirical determinations with appro- 
priate priors derived from theoretical models or observational 
constraints. Several authors began to produce 3D reddening 
maps using optical and/or near-infrared photometry and spec- 
troscopy from large surveys adopting different assumptions on 


dust and stellar distributions 
2010). Among these, we chose the one recently presented in 
Potent rjr MM of its accurate description of the 
local interstellar medium and its overall spatial coverage which 
allowed us to determine the reddening for a large fraction of PIC 
stars in a homogeneous way. The 3D reddening map of 
is described in Appendix[E] The median red- 
dening of asPIC1.1 stars inside the reddening map is equal to 
E(B — V)=0.04 and the median uncertainty is E(B — V)=0.02. 
For the M sample, 99.8% of the stars are contained in the redden- 
ing map and, for the FGK sample, this figure is 81.8%. For the 
stars falling outside the map, reddening is calculated up to the 
edge of the map, and then a correction is added, as described in 
Sect.[E.1] Figure[11 presents the E(B-V) versus distance diagram 
for all stars in asPIC1.1. 


5. Stellar parameters 


This section describes the algorithms we used to estimate effec- 
tive temperatures, radii, and masses. 


8 Very recently |Lallement et al.|(2019) published an updated version 


of the map, which is not included in this version of the catalogue but 
will be considered in the next version. 
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Fig. 11. Reddening E(B-V) vs distance for all stars in asPIC1.1. 


5.1. Intrinsic colour-effective temperature relation 


We first determined the relationship between the intrinsic colour 
in the Gaia bands and the effective temperature. For FGK spec- 
tral types, we used the stellar sample from {Casagrande et al.| 
(2010). Effective temperatures reported in this catalogue are 
calculated using the infrared flux method (IRFM). We cross- 
matched this sample with Gaia DR2 and considered only stars 
closer than 40 pc from the Sun. (2010) derived 
the IRFM effective temperatures for all stars with HIPPARCOS 
parallaxes and distances € 70 pc from the Sun assuming they are 
unaffected by reddening. To be consistent with this choice we 
assumed that the sample of calibration stars we used (all located 
within 40 pc) have null reddening. The astrometric and photo- 
metric quality conditions expressed by Eq. (1) and Eq. (2) in 
were verified for all the selected starg?| Our 
final list comprises 110 stars out of the 423 stars of 
etal. ](2010). 


For M-dwarfs, we considered the sample of[Mann et al.|(2015). 
All the 179 (see Sect. [3.2 stars in the sample are located within 


40 pc from the Sun. We also checked that Eq. (1) and Eq. (2) of 
were satisfied, selecting in this way a total of 
171 stars. Their effective temperatures were estimated from op- 
tical and near-infrared spectra. Finally, for hot stars (Teg >6510 
K), we considered the temperature and colour estimates reported 


in|Pecaut & Mamajek| (2013). 


As shown in Figure from the resulting sample of 281 
stars we interpolated a single fifth-order polynomial relating the 
observed Gaia colour of the calibration stars and their effective 
temperature. The coefficients of the fit are reported below, and 
the RMS of the residuals of the fit is equal to 65 K 


? With the exception of one source, the bright spectroscopic binary 
HD112758, for which the astrometric condition expressed by Eq. (1) of 
was not satisfied, but that was kept in the sample 
because it is known to host a low-stellar-mass companion 


Quirrenbach|201 1). 
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Fig. 12. Relation between effective temperature Ter and intrinsic colour 
(Ggp — Ggp)o. The continuous black line denotes our best-fit model 
(Eq. 4). while the dashed vertical lines denote the limits of validity of 
the relation. Coloured dots represent the samples of |Casagrande et al. 


(2010), (2015), and|Pecaut & Mamajek|(2013) discussed in 
the text in Sect. 


T(K) = 9453.14 — 6859.40 (GBP = Grp)o+ 
+ 3542.16 (Gsp — Grp)2+ 
— 1053.09 (Gpp — Grp)? + 165.635 (Gpp — Grp)4+ 


— 10.5672 (Gap — Grr), (4) 


where the relation is valid for 0.5 <(Ggp — Grp) < 5. Hereafter, 
we use Eq. as the relation between the intrinsic colour and 
the effective temperature. Using Eq. (4) we obtained a first guess 
of the effective temperature of all asPIC stars. These stars are 
generally located further, at larger distances than the sample of 
calibration stars, and therefore reddening cannot be neglected. 
As reported in Section |4| we estimate reddening using the 
reddening map Men a (2018). The customised algo- 
rithm described in Section. |E|furnishes the interpolated redden- 
ing E(B — V) for any 3D star position. As input, it requires the 
Galactic longitude / and the Galactic latitude b which are taken 
directly from the input catalogue along with the distance from 


the|Bailer-Jones et al.|(2018) catalogue. 


5.2. Conversion from E(B — V) to E(Ggp — Grp) and 
determination of extinction Ag 


The E(B — V) reddening is converted to E(Gpp — Grp) and to the 
extinction Ag in the Gaia G-band as follows: 


Ag = Rg E(B - V), (5) 
App = Rpp E(B —- V), (6) 
App = Rgp E(B — V), (7) 


where the Gaia extinction coefficients depend on the effective 
temperature as follows 


Rg = — 0.5335 + 12.9373 (T4) — 13.9514 (T4) + 


— 13.8012 (T4)? + 40.9902 (T4)* — 23.6648 (T4)°, (8) 
Rpp = — 2.4689 + 59.5802 (T4) — 253.9922 (T4 + 

+ 526.5333 (T4)? — 523.9970 (T4 + 

+ 201.2829 (T4), (9) 
Rgp =0.0407 + 9.8825 (T4) — 16.8207 (T4)?+ 

+ 9.1283 (T4)? + 2.5051 (T4)! — 2.4158 (T4)°, (10) 
and 
T4 = 1074 (Teg). (11) 


Such relations were derived using Padova stellar models 
[san et al.|[2012), assuming solar metallicity, solar composition, 
and temperatures within the interval 2500 K < Ter < 7200 K. In 
particular, we modelled the ratio of the absorption in a given pho- 
tometric band to the reddening E(B — V) considering an unred- 
dened 1 Gyr isochrone, and the same isochrone reddened us- 
ing Ay=0.5 mag, where reddening was applied on a star-by-star 
basis. We then interpolated apparent magnitudes versus effec- 
tive temperature on the same temperature grid for both the red- 
dened and unreddened isochrones, and calculated the absorption 
on each Gaia band as a function of temperature. In Ps we 
show our interpolating relations presented in Eqs. (8). (9), and 
(10) along with the results obtained using the bolometric correc- 
tion program of{Casagrande & VandenBerg|po1s]"] In particu- 
lar, we used this program to calculate the bolometric corrections 
for unreddened and reddened stars interpolating across the grid 
of log g and effective temperatures of our adopted isochrone (as- 
suming solar metallicity and solar composition), and from the 
bolometric corrections we calculated the ratios between the ex- 
tinction in each band and the reddening E(B-V) as a function of 
the effective temperature| | The results we obtained from our 
procedure agree well with those of 
in general. Some small differences are visible, partic- 
ularly for late-type stars and for the G and Ggp bands which 
may be explained by differences in the underlying model at- 
mospheres, because the Padova group colours and magnitudes 
are based on the libraries of stellar spectra of 


(2003), whereas |Casagrande & VandenBerg| (2018) used 
the MARCS library of theoretical stellar fluxes (Gustafsson et al. 


22008}. 


5.3. Correction for reddening and extinction 


The observed Ggp— Grp colour and G magnitude were corrected 
for reddening and extinction to obtain the intrinsic colour (Gpp — 
Grp)o and the intrinsic magnitude Go : 


(Gap — Gnp)o = (Gpp — Grp) - E(Gpp — Grp), (12) 


Go =G-Ag. (13) 


10 https://github.com/casaluca/bolometric-corrections 


11 Ay |. (BCay-9-BCay) = _ 
EG sy x3.1, where y = Ggp,G,Grp and Av=0.5. 


Article number, page 11 of 23 


A&A proofs: manuscript no. PicPaperCorr 


36 ————————Ó-— 


| _ Casagrande 


A,/E(B- V) 


— — Casagrande & VandenBerg (2018) | 


& VandenBerg (2018) 


Casagrande & VandenBerg (2018) 


| LEE ride n E EE ee Ae a ea rral 


ia V i4. TATE 
4000 5000 6000 7000 


Tore (K) 


Li 1 
3000 


6000 7000 4000 5000 6000 7000 


Tog (K) 


Fig. 13. Relations between Ag/E(B — V) (left), Ac,, /E(B — V) (middle), Ac, / E(B — V) (right) and effective temperature Ter derived with the 
procedure described in the text (Sect. B2). The red lines represent our interpolating relations and the blue lines models from |Casagrande & 


VandenBerg (2018). 


5.4. Effective temperature: the iteration cycle 


We recalculated the effective temperature using Eq. and the 
new intrinsic colour (Gpp — Ggp)o in Eq. (12). We then used the 
new effective temperature to estimate the reddening coefficients 
(Section[5.2) and then a new estimate of the intrinsic colour (Sec- 
tion[5.3). We then used this to estimate the effective temperature 
again. A subsequent iteration lead to effective temperatures that 
differ by less than 10 K from last estimate; therefore two itera- 
tions were sufficient. 


5.5. Bolometric correction in the G—band, BC 


The bolometric correction of the G-band is obtained from the 


relations given in (2018): 


BCe = 6.0 x 10? + 6.731 x 10? AT og — 6.647 x 1095 AT Z+ 


+ 2.859 x 107! ATS. — 7.197 x 107 ATA, (14) 
for a temperature interval 4000 K < Ter < 8000 K and 
BCg = 1.749 + 1.977 x 10? AT og + 3.737 x 1077 AT2,- 

— 8.966 x 101! AT3, — 4.183 x 10 AT fp, (15) 


for a temperature interval 3300 K < Teg < 4000 K and where 
AT. = Tet = Tefo and Tefo = 5772K. 


5.6. Determination of the absolute magnitude Mg o and of 
the absolute luminosity L 


The intrinsic absolute magnitude Mg and the luminosity L were 
calculated as follows: 


Mao = Go -5 logio d+ 5. (16) 


L = 10794 Ma.o+BCa-Mroro) 
Lo i 


where d is the distance of the star in parsecs from Bailer-Jones| 
let al.|(2018) and Mpox, = 4.74 and Lo = 3.828 x 107 W. 
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(17) 


5.7. Determination of the stellar radius 


We then calculated the stellar radius R from the Stefan- 


Boltzmann law: 


(re) 


The presence of absorption and reddening influences the deter- 
mination of the stellar radius because it affects both the estima- 
tion of the effective temperature and of the absolute luminosity. 
Therefore it is important to apply dereddening procedures like 
those we describe above. Reddening and absorption have two 
opposite and partially compensating effects on the stellar radius. 
A more detailed discussion about this point is provided in Ap- 


pendix 


Tef 
Tefo 


L 
Lo 


(18) 


5.8. Determination of the stellar mass 


The stellar mass was determined using the following relation 


from|Moya et al.|(2018): 


M L 
toe 2] = -0.119+2.14x 10? T. +0.1837 log £). (19) 
Mo Lo 


The range of validity is 4780 K < Teg < 10990 K and -0.717 < 
log (L/Lo) < 2.01. This relationship is valid for dwarfs and sub- 
giant stars. 


5.9. Impact of metallicity on the determination of stellar 
parameters 


The procedure used to determine stellar parameters described in 
Sect.[5]does not take into account metallicity because we do not 
have metallicity measurements for most of our stars. In this sec- 
tion, we estimate the impact that metallicity has on the determi- 
nation of stellar parameters. We retrieved a grid of isochrones 
with age between 7«log(Age/yr)«10 in steps of 0.2 dex and 
metallicity between -1«[M/H]«0.3 in steps of 0.1 dex from the 
Padova stellar database. For each model in the grid, we deter- 
mined the stellar parameters with the procedure described in 
Sect. [5] and compared them with the theoretical parameters re- 
ported in the isochrones for six different intervals of metallicity. 
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Fig. 14. Relative differences between stellar parameters estimated using asPIC1.1 and theoretical parameters from stellar isochrones in different 
metallicity intervals. The rightmost labels of each panel denote the centres of the [M/H] metallicity intervals having a width of +0.1 dex as also 


listed in Table] 


Table 3. Median relative differences (as a percentage) and standard de- 
viations between asPIC1.1 estimated stellar parameters and theoretical 
stellar parameters in different metallicity intervals. 


[M/H] aT (%) RAO My) 
-0.9<[M/H]<-0.7 . 0.02x0.33 -2.5+0.4 1749 
-0.7<[M/H]<-0.5 -0.04+0.30 -1.840.3 1447 
-0.5<[M/H]<-0.3 -0.06+0.33 -1.10.4 11+6 
-0.3<[M/H]<-0.1 -0.21+0.38 -0.10.5 6+5 
-0.1<[M/H]<0.1 -0.41+0.36 1.0+0.5 0.5+3.9 
0.1<[M/H]<0.3 -0.66+0.36 2.6+0.8 -4.143.5 


In particular, we calculated the relative differences ( ar, AS, aM 


between our parameters (Te aspici.1, Radiusaspici.1, Massaspici.i1) 
and the theoretical parameters (Tef Iso, Radius;so, Massrso) 


AT — Tefraspici.1 — Teft.so 


(20) 
T Tef, Iso 
AR _ Radiusaspici.1 — Radius;so (21) 
R Radiusiso 
AM m MasSasPIct.1 — Massiso (22) 


M Massjso 


The results are presented in Fig. [14] and in Table 3] We find 
that the effective temperature has a weak dependence on metal- 
licity with relative differences ranging from -0.6% to 0.02%. 
The radius presents relative differences of between -2.5% and 
2.5% while the mass shows the largest relative differences rang- 
ing from -4% to 17%. We note that the mass is also the most 
uncertain among the theoretical parameters associated to the 
isochrones. Both the temperature and the mass relative differ- 
ences are negatively correlated with the metallicity while the 
radius is positively correlated. We also calculated the standard 
deviations of the relative differences in each metallicity interval. 
These differences range between 0.3% and 0.4% for the temper- 
ature, between 0.3% and 0.8% for the radius, and between 4% 
and 9% for the mass. 


5.10. Errors estimate 


We determined the errors on stellar parameters using Monte 
Carlo simulations, perturbing all observing quantities accord- 
ing to their associated errors, and assuming a Gaussian pertur- 
bation. Beyond colours, magnitudes, distances, and reddening, 
we also perturbed the bolometric corrections adopting the errors 


reported by (2018). The error on the reddening was 
taken from the reddening map of|Lallement et al.|(2018) for stars 


falling inside its boundaries. For stars falling outside this map, 
the error on the reddening was assumed equal to twice the er- 
ror on the median reddening of all stars inside the map, which is 
equal to a E(B — V)=0.04 (see also Sect. [E). The effective tem- 
perature deduced by the colour-effective temperature relation 
was perturbed considering a conservative error of 200 K on its 
determination. We performed 100 simulations for each star and 
then calculated the errors on the effective temperature, mass, and 
radius as the half interval between the 16" and the 84" percentile 
of the cumulative distribution of the simulated values. The me- 
dian value of the effective temperature error is 227 K (4%), and 
this value is 0.1 Re (7%) for the radius, and 0.1 Mo (8%) for the 
mass. 


5.11. Stellar parameter distributions 


In Fig. [I5] we present the distributions of effective temperature, 
radius, and mass for the FGK sample. We note that stellar pa- 
rameters have not yet been calculated for the M dwarfs sample. 
They will be included in the next version of the catalogue. 


6. Comparisons 


We compared stellar parameters and reddening in asPIC1.1 with 
the values reported in the TIC (Stassun et al.|2019). We used 
the Candidate Target List v8.01 (CTL v8.01) and cross-matched 
it with asPIC1.1 using the Gaia source ID reported in the CTL 
catalogue. We considered only those stars for which all stellar 
parameters and the reddening are defined in both catalogues. 
We identified 2 022 913 stars in common between the two cat- 
alogues. The median differences and the standard deviations of 
the differences between asPIC1.1 and TIC temperatures, radii, 
masses, and reddenings are: AT.g = (—100 + 300) K, AR, = 
(0.05 + 0.07) Ro, AM, = (0.05 + 0.17) Mo, and AE(B - V) = 
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Fig. 15. Distributions of effective temperature (left panel), radius (middle panel), and mass (right panel) for the FGK sample. 


Catalogue AT ot AR, AM, Alog g AE(B — V) 
(K) (Ro) (Mo) (dex) (mag) 
TIC (CTLyg) -100+300 0.05+0.07 0.05 x0.17 - —0.002 + 0.055 
GALAH DR2 4100 - - 0.10.2 - 
HARPS -70+90 - - -0.03+0.15 - 
APOKASC (SDSS) -130+90 0.10.2 0.10.1 -0.02+0.06 - 
APOKASC (ASPCAP) 90+130 0.1+0.1 0.1+0.1 -0.01+0.06 - 


Table 4. Median differences and standard deviation of the differences between asPIC1.1 and the stellar parameters and interstellar reddening values 


from other catalogues. 


(—0.002 + 0.055). We note that the distribution of stellar masses 
appears markedly asymmetric with asPIC having a tendency to 
overestimate the mass with respect to the TIC. This fact may 
arise from the different choices adopted for the calibration of 
stellar masses in the two catalogues. In our case, we adopted the 


calibration from (2018) which accounts both for the 


temperature and absolute luminosity of the stars, while Stassun| 
and [Stassun et al.|(2019) adopt a calibration based 
only on the temperatures. The differences between the asPIC1.1 
and TICv8 estimated parameters are shown in the histograms in 
Fig. [16] 

We then considered the Galactic Archaeology with HER- 
MES (GALAH) survey second data release which contains 
342 682 spectroscopically observed stars (Buder et al.[2018). We 
cross-matched GALAH with asPIC1.1 using the Gaia source ID 
reported in the catalogue. We found 64 061 matched sources. We 
compared the effective temperatures and the log g in GALAH 
with the asPIC1.1 values. In asPIC1.1 we report stellar masses 
(M) and radii (R). The log g is estimated using log g = log M 
- 2 log R + 4.4374 (Smalley|2005). The median differences be- 
tween the asPIC1.1 and the GALAH effective temperatures and 
gravities are: AT,g = (4 + 100) K and Alog g = (0.1 + 0.2) dex, 
as shown in Fig. 

We also considered the sample of 1111 FGK dwarf stars 
from the HARPS GTO program, homogeneously analysed in 
(2012). We matched the catalogue with as- 
PICI.1 using angular distances (accounting for proper motions) 
and found 1102 common sources. The median differences be- 
tween the asPIC1.1 and the HARPS effective temperatures and 
gravities are: AT. = (—70 + 90) K and Alog g = (—0.03 + 0.15) 
dex (Fig. [18). 

Finally, we compared our parameters with the ones reported 
in the first APOKASC catalogue of spectroscopic and asteroseis- 


mic data for dwarfs and subgiants (Serenelli et al.|2017). The 
cross-match with asPIC1.1 yielded 385 common stars for which 


Article number, page 14 of 23 


all parameters were defined in both catalogues. The analysis was 
performed both for the SDSS temperature scale and for the spec- 
troscopic temperature scale adopted in the APOGEE Stellar Pa- 
rameters and Chemical Abundances pipeline (ASPCAP). In the 
first case, the median differences and the standard deviations 
of the differences between asPIC1.1 and APOKASC tempera- 
tures, radii, masses, and gravities are: AT.g = (—130 + 90) K, 
AR, = (0.1 + 0.2) Ro, AM, = (0.1 + 0.1) Me, and Alogg = 
(—0.02 + 0.06), while in the second case: ATeg = (90 + 130) K, 
AR, = (0.1 + 0.1)Ro, AM, = (0.1 + 0.1) Me, and Alogg = 
(—0.01 + 0.06) (see also Fig.[19). 


Table [4] summarises the results of the comparisons between 
asPIC1.1 and the stellar catalogues considered in this section and 
we conclude that the average difference between our effective 
temperatures and those of the other catalogues is -40 K, the av- 
erage radius difference is 0.08 Ro, and the average mass differ- 
ence is 0.08 Mo. Considering the internal errors we estimated in 
Sect. [5] and adding them in quadrature to the estimated external 
errors coming from the comparisons from other catalogues (see 
above), we conclude that our overall (internal+external) uncer- 
tainties on the stellar parameters determination is 230 K (4%) 
for the effective temperatures, 0.13 Ro (9%) for the stellar radii, 
and 0.13 Mo (11%) for the stellar masses. 


7. Special target list 


With the asPIC 1.1 catalogue we also release a special list of 
objects that consists of all currently known confirmed planet 
hosts included in the asPIC1.1. This list has been constructed 
using the Exo-MerCat tool (Alei et al.|[2020), which accesses 
four main public catalogues of exoplanets, and obtain an up- 
dated list of the known planets along with their physical and 
orbital parameters. In particular, this tool considers the Extra- 
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Fig. 16. Distributions of effective temperature, stellar radius, mass, and interstellar reddening differences between asPIC1.1 and CTL. 


solar Planets Encyclopaedid"] the NASA Exoplanet Archivd] 
the Exoplanets Orbit Databasq"] and the Open Exoplanet cat- 
alogud"] each one employing its own cataloguing system. For 
this reason, before obtaining a uniform merging of the objects of 
interest, Exo-MerCat performs a proper discrimination among 
aliases and a standardisation of the different entries provided by 
the catalogues. The list of the known exoplanets is matched with 
the asPIC1.1 catalogue. 


http://exoplanet.eu/ 


ttps://exoplanetarchive.ipac.caltech.edu/ 


ttp://exoplanets.org 


ttp://www.openexoplanetcatalogue.com/ 


The special targets list is a living catalogue: it will be contin- 
uously updated at any new release of asPIC, because the num- 
ber of newly discovered planets will constantly increase over the 
years, not only before but also during the PLATO mission life- 
time. 


8. Conclusions 


In this paper, we present the asPIC1.1 catalogue, a public all-sky 
catalogue of dwarf and subgiant stars of interest for the PLATO 
survey, based on the Gaia DR2 data release. The asPIC cata- 
logue will be fundamental to identifying the best fields for the 
PLATO space mission and the most promising targets, analysing 
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Fig. 17. Distributions of effective temperatures, and stellar gravity differences between asPIC1.1 and GALAH/DR2. 
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Fig. 18. Distributions of effective temperatures, and stellar gravity differences between asPIC1.1 and the HARPS sample. 


the instrumental performances as well as planning and optimis- 
ing ground-based follow-up studies. This catalogue also repre- 
sents a valuable resource for the construction of stellar samples 
optimised for the search of transiting planets. 


The catalogue includes a total of 2 675 539 stars among 
which 2 378 177 are FGK dwarfs and subgiants and 297 362 
are M-dwarfs. The median distance of FGK stars in our sample 
is 428 pc and that for M dwarfs is 146 pc. 


Article number, page 16 of 23 


We also show that our selection criteria do not bias the sta- 
tistical distribution of metallicities, and we analyse the impact 
that metallicity has on the derivation of stellar parameters. We 
derived the reddening of our targets and developed an algo- 
rithm to infer stellar fundamental parameters from astrometric 
and photometric measurements. We show that the overall (inter- 
nal+external) uncertainties on the stellar parameters determined 
by our analysis are ~230 K (4%) for the effective temperatures, 
~0.1 Re (9%) for the stellar radii, and «0.1 Mo (11%) for the 
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Fig. 19. Distributions of effective temperatures, stellar radii, masses and gravities differences between asPIC1.1 and APOKASC sample and two 


different temperature scales (SDSS, black and ASPCAP, red) . 


stellar masses. We also release a special target list containing all 
known planet hosts cross-matched with our catalogue. 
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Appendix A: Johnson and Gaia DR2 colour 
transformation 


The definition of PLATO samples requires knowledge of the vi- 
sual apparent magnitude V. The V magnitude we used to select 
the targets comes from a transformation of Gaia DR2 magni- 
tudes and colours to the Johnson V band. We decided to adopt 
Gaia DR2 colours and magnitudes as the basis of our calcula- 
tion of the V band magnitude because Gaia offers a homoge- 
neous, precise, deep, multi-band, all sky photometry.[Evans et al.] 
(2018) give colour-colour relations between several commonly 
used passbands and Gaia DR2 colours. However, the colour 
range within which these relationships are applicable does not 
include all the spectral types of interest for PLATO, in particular 
the M-dwarfs. In order to extend the|Evans et al.|(2018) relation- 
ship to the entire M-type regime, we used a sample of main se- 
quence stars covering a broad range of spectral types for which 
accurate and precise photometry in the Johnson V is given. To 
this aim, we used photometry of standard stars given in the cat- 
alogues of and [Landolt] (2009). In addition, for 
M-dwarfs we used the RECONS catalogue (Henry et al./2018) 
which provides an accurate census of a volume-limited sample 
of neighbour stars and is a valuable resource in particular for M 
dwarfs which are poorly represented in the other catalogues. The 
Stetson, Landolt, and RECONS catalogues were cross-matched 
with Gaia DR2. We then corrected each star for reddening using 
the procedure described in Sect. 4| and isolated main sequence 
stars within 200 pc from the Sun with error on the Gaia G-band 
photometry og < 0.003, reddening E(B-V)«0.05, and satisfy- 
ing the astrometric and photometric quality criteria defined by 
Eq. (1) and Eq. (2) of|Arenou et al.| (2018) ). In total we selected 
267 stars from|Stetson (2000), 98 sis Toca emat 09) and 
144 stars from (2018). We interpolated the (G-V)o 
versus (Gpp-Grp)o colours with a fifth-order polynomial relation 
and the best-fit coefficients are reported in Table D] of the main 
text. Figure|[A. I|shows a graphical representation of the relation- 
ship, overlaying the best-fit model (solid red line) and the stellar 
samples that have been used (coloured dots), the 
relation (solid black line), and the (2019) 
calibration (dashed black line) used in the TICv8, and based on 
PHOENIX model atmospheres (provided by K. Stassun, private 
communication). The standard deviation of the residuals from 
our best-fit model is oj =0.02745 mag. The relation is valid 
within 0.5<(Ggp-Grp)o «5. We note that this relation is applied 
to the intrinsic (i.e., reddening-free) magnitudes, and then ab- 
sorption is added (see Sect.H) to obtain the apparent magnitude. 


Appendix B: Selection of FGK dwarfs and 
subgiants using Johnson B,V photometry. 


For completeness of our discussion we report here an alter- 
native selection of FGK dwarf and subgiant stars used during 
the construction of the first version of the catalogue (named 
asPIC1.0). The procedure is similar to the one described in 
Sect. 3| but is based on the Johnson photometric system. The 
selection criteria are then defined in the absolute and intrinsic 
Myo versus (B — V)o colour magnitude diagram and are given by 
the following inequations where we further distinguish dwarfs 
from subgiants: 
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Fig. A.1. Relationship between the intrinsic (G-V)ọ and (Ggp-Gpp)o 
colours calibrated in this work (solid red line). We also show the corre- 


sponding [Evans et al.] (2018) (solid black line) and|Stassun et al.] (2019) 


(dashed black line) relationships. 


Dwarfs: 


0.42 < (B — V)o < 1.38 
Myo = 5(B — V)o + 0.4 
Mvo «5(B- V); - 3.5 
V «13 


Subgiants: 


0.42 < (B — V) € 0.8 
Myo < 5(B — Vo + 0.4 
Myo > 5(B- V)o — 2] 
V «13 


OR 


0.8 < (B- V); € 1.0 
Mvo «4.5 

Myo > 5(B- V)o — 2 
V < 13. 


The selection is depicted in Fig. One of the limitations of 
this selection is that it cannot be easily extended in the M-dwarf 
domain because of saturation effects of the (B — V) colour. With 
the release of Gaia DR2 photometry we decided to work directly 
in the Gaia photometric system, as described in Sect. |3| Never- 
theless the relationships reported here can be useful to isolate 
FGK dwarfs and subgiants using Johnson B and V photometry. 
The TPR and FPR of this selection are respectively 10096 and 
45%, and therefore such a selection has the same effectiveness 
as that presented in Sect. 3.3] in recovering true positives, but 
implies a larger rate of false positives. 
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Table B.1. Description of the columns of the asPIC1.1 catalogue. 


Col. number Col. name Units Description 

1 PICidDR1 - PLATO identifier for asPIC1.1 

2 PICnameDR1 - PLATO identifier for asPIC1.1 with version number 

3 sourceld - Gaia DR2 sourceld 

4 RAdeg deg Barycentric Right Ascension (o) of the source in ICRS at the reference epoch 

5 eRAdeg mas Standard error Ce» = 0,cosó of the right ascension in ICRS at the reference epoch 

6 DEdeg deg Barycentric Declination (6) of the source in ICRS at the reference epoch 

7 eDEdeg mas Standard error c; of the declination in ICRS at the reference epoch 

8 Plx mas Parallax 

9 ePlx mas Standard error on parallax 

10 pmRA mas yr! Proper motion in Right Ascension direction He» = Hacosd in ICRS at the reference epoch 
11 epmRA mas yr Standard error on proper motion in right ascension direction 

12 pmDE mas yr! Proper motion in Declination direction jis in ICRS at the reference epoch 

13 epmDE mas yr! Standard error on proper motion in declination direction 

14 PM mas yr? Total proper motion calculated as PM= ./pmRA2 + pmDE? 

15 ePM mas yr! Error on the total proper motion 

16 Epoch yr Reference epoch to which the astrometric source parameters are referred 

17 Gflux e s! Mean flux in the G band 

18 eGflux e s! Error on G band mean flux 

19 Gmag mag Mean magnitude in the G band 

20 eGmag mag Error on mean magnitude in the G band 

21 BPflux e s Mean flux in the integrated BP band 

22 eBPflux e7 s7! Error on mean flux in the integrated BP band 

23 BPmag mag Mean magnitude in the integrated BP band 

24 eBPmag mag Error in the mean magnitude in the integrated BP band 

25 RPflux e s! Mean flux in the integrated RP band 

26 eRPflux e7 s7! Error on mean flux in the integrated BP band 

27 RPmag mag Mean magnitude in the integrated RP band 

28 eRPmag mag Error in the mean magnitude in the integrated RP band 

29 BPRPexcess - BP/RP excess factor 

30 BPRP mag Gpp — Gpp colour 

31 BPG mag Gpp — G colour 

32 GRP mag G — Gpp colour 

33 GLON deg Galactic longitude of the object at the reference epoch 

34 GLAT deg Galactic latitude of the object at the reference epoch 

35 ELON deg Ecliptic longitude of the object at the reference epoch 

36 ELAT deg Ecliptic latitude of the object at the reference epoch 

37 rest pe The estimated distance 

38 rlo pe Lower bound of the confidence interval of the estimated distance (Bailer-Jones et al.|2018) 
39 rhi pc Higher bound of the confidence interval of the estimated distance (Bailer-Jones et al./20T8) 
40 erest pc Error on the estimated distance calculated as restError= 0.5(rhi — rlo 

41 AG mag Estimate of the extinction Ag in the G band 

42 eAG mag Uncertainty on Ag estimate 

43 EBPRP mag Estimate of reddening E(Ggp - Grp) 

44 eEBPRP mag Uncertainty on reddening E(Ggp - Grp) 

45 EBV mag Estimate of reddening E(B-V) 

46 eEBV mag Uncertainty on E(B-V) reddening 

47 extStatus - extStatus value. If extStatus=1 the star is outside the reddening map of[Lallement et al.|(2018} 
48 BPRPO mag Dereddened (Gpp - Grp) colour 

49 eBPRPO mag Error on dereddened (Gpp - Grp) colour (BPRPO) 

50 gaiaVO mag Dereddened apparent visual magnitude V in the Johnson-Cousin system 

51 egaiaVO mag Error on intrinsic visual magnitude (gaiaVO) 

52 BJgaiaMVO mag Absolute intrinsic visual magnitude obtained as: BJgaiaMV0-gaiaVO-5log,, rest + 5 
53 eBJgaiaMVO mag Error on absolute intrinsic visual magnitude (BJgaiaMVO) 

54 gaiaV mag Visual apparent magnitude in the Johnson-Cousin system obtained as: gaiaV=gaiaV0+3.1 EBV 
55 egaiaV mag Error on apparent visual magnitude V (gaiaV) 

56 BJgaiaMGO mag Absolute intrinsic G magnitude obtained as: BJgaiaMG0-Gmag-5log,, rest + 5 - AG 
57 eBJgaiaMGO mag Error on intrinsic absolute G magnitude (BJgaiaMGO0) 

58 Teff K Estimate of stellar effective temperature (Ter) 

59 eTeff K Uncertainty on T. estimate 

60 Radius Ro Estimate of stellar radius 

61 eRadius Ro Uncertainty on stellar radius estimate 

62 Mass Mo Estimate of stellar mass 

63 eMass Mo Uncertainty on mass estimate 

64 sourceFlag - Bitmask with the following meaning: 


sourceFlag=1 - FGK star 
sourceFlag=2 - M star 
sourceFlag=4 - Planet Host 
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Table B.1. - Continued. 


Col. number Col. name Units Description 

65 contGaiaMag60 mag Total G magnitude of the Gaia contaminants within 60 arcsec 
66 contNumber60 - Total number of Gaia contaminants within 60 arcsec 

67 contGaiaMag45 mag Total G magnitude of the Gaia contaminants within 45 arcsec 
68 contNumber45 - Total number of Gaia contaminants within 45 arcsec 

69 contGaiaMag30 mag Total G magnitude of the Gaia contaminants within 30 arcsec 
70 contNumber30 - Total number of Gaia contaminants within 30 arcsec 


Mvo 


10 


05 06 07 08 09 10 11 12 13 
(B — V)o 


Fig. B.1. Selection of FGK dwarf and subgiant stars in the Johnson B, 
V absolute intrinsic colour-magnitude diagram. 


Appendix C: The catalogue 
Appendix C.1: Naming convention 


The asPIC filename (asPIC.n1.n2) contains two dot-separated in- 
teger numbers which are used to distinguish the different ver- 
sions of the catalogue. The first number (n1) indicates an asPIC 
major release and is linked to the Gaia data release used to build 
the catalogue. When a new Gaia release triggers an asPIC major 
update, nl is increased by one. The second number (n2) indi- 
cates a substantial update, which is not major and thus does not 
require changing n1. Examples of substantial updates are, for ex- 
ample, a change in (a) the number of objects, (b) the number of 
columns, (c) the format of one or more columns, or (d) the con- 
tent of one or more columns. When the first number n1 changes, 
n2 is set to zero. 


Appendix C.2: Content 


In Table [B.1| we report the asPIC1.1 column names and a brief 
description of their meaning. 


Appendix D: Implementation procedure 


In this Appendix we describe the implementation procedure 
adopted to build the asPIC1.1 catalogue. The implementation 
of PLATO asPIC1.1 is done using the PLATO Data Processing 
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Centre - ASI (PDPC-A) infrastructure dedicated to the Gaia and 
PLATO missions. 


— The procedure starts using the Gaia DR2 catalogue main 
table (gaiaSource) and the Gaia estimated distances table 
from (2018). A asPIC main table (as- 
PIC1.1) is created and all sources from gaiaSource where: 
photGMeanMag x 15.5 AND bpRp is not null AND rest is 
not null are added. The asPIC1.1 table includes 52 557 711 
sources. 

— The asPIC1.1 table is updated and the following columns are 
added: pmtotal, pmtotalError, photGmeanMagError, 
photBpMeanMagError, photRpMeanMagError, 
restError. 

— We created an input file to be fed into the stellar parame- 
ters’ code (Sect. [5). The input file includes all sources in 
the asPIC1.1 table. The code is run in parallel. The results 
of the stellar parameters’ code (47 966 591 sources) are 
then fed back into the database in a dedicated table. The 
stellar parameters’ code output columns are: ag, agError, 
ebprp, ebprpError, ebv, ebvError, extStatus, teff, 
teffError, radius, radiusError, mass, massError and 
dwsgFlag. The flag dwsgFlag discriminates between FGK 
(dwsgFlag=1), M (dwsgFlag=2) stars and stars that are nei- 
ther FGK or M (dwsgFlag=0). The main table (asPIC1.1) is 
updated to include these results. 

— The main table is further updated to include bpRp9, 
bpRpOError, BJgaiaMGO0, BJgaiaMGOError, PICidDR1, 
gaiaV6, gaiaV8Error, gaiaV, gaiaVError and gaiaMV9, 
gaiaMVOError. Only 47 288 759 sources out of 47 966 591 
are updated, since a few sources are left out because they do 
not fall within the Gaia-Johnson calibration limits. 

— The definition of the samples in the asPIC110 table is up- 
dated according to the final selection criteria: 

1. FGK: (dwsgFlag-1 AND  gaiaV«-13.0 AND 
BpRp0>0.56) which outputs 2 378 177 sources; 

2.M: (dwsgFlag=2 AND gaiaV<=16.0 AND 
(BJgaiaMGO >  2.334"(BpRp8) + 2.259) AND 
rest<600.0) which outputs 297 362 sources. 

— Photometric contaminants are searched for and counted in 
the Gaia DR2 catalogue for all the candidate targets listed in 
asPIC using a cone search with 30, 45, and 60 arcsec radius 
respectively. To compute the total G magnitude we summed 
the entire flux of each contaminant. As the contaminants 
may fall close to the circle border, the computed total 
magnitude is over-estimated. This information was used to 
compute the columns contGaiaMag30, contNumber390, 
contGaiaMag45, contNumber45,contGaiaMag60 
contNumber69. 


and 


Appendix E: The 3D reddening map 
The 3D reddening map presented in |Lallement et al.| (2018) 


is produced using a hierarchical approach merging individual 
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photometric colour excess measurements of stars with colour 
excesses estimated from diffuse interstellar bands absorption. 
These measurements are coupled with a non-homogeneous, 
large-scale prior reddening distribution deduced from massive 
photometric surveys. The 3D map eA 
was retrieved from the dedicated websitq!'| The map is extended 
in a region between —2kpc and 2kpc from the Sun position 
in the Heliocentric Galactic X and Y coordinates and between 
—300 pc and 300 pc in the Galactic Z coordinate. It has a reso- 
lution of 5 pc? and units are mag - pc^!. The map is particularly 
useful at small distances since it reveals how local cavities and 
cloud complexes are spatially distributed around the Sun. 


We used the map of|Lallement et al.|(2018) in conjunction with 
Bailer-Jones et al.| (2018) distances from Gaia DR2 to esti- 


mate the amount of interstellar reddening to our target stars. We 
evaluated the reddening E(B — V); on a sample of N, equally 
spaced, nodal points along each target direction by considering 
the weighted average of the 27 surrounding voxels reddening (a 
3 x 3 x 3 voxels cube around the nodal point). The reddening 
was weighted by the distance of the nodal point from each one 
of the neighbouring voxels' centers while the distance between 
two adjacent nodal points was 2 pc. We therefore summed up the 
reddening along the boresight to derive the integrated reddening 
to the target star 


i=N 


dE(B — V) 
E(B-V)- 5| — p AR (pe), (E.1) 
il 
where 
AR = 2pc 
and 
dE(B — V)i 
dR 


is the reddening per unit lengtH['3] in the i-th voxel. The median 
reddening of asPIC1.1 stars inside the reddening map is equal to 
E(B — V)=0.04 and the median uncertainty is c E(B — V)=0.02. 
For the M sample 99.846 of the stars are contained in the redden- 
ing map and, for the FGK sample, 81.8%. For the stars falling 
outside the map, reddening is calculated up to the edge of the 
map, and then a correction is added as described in Sect. [E.T] 


Appendix E.1: Reddening map extension 


As discussed in the previous section, a small fraction of asPIC 


stars falls outside the limits of the|Lallement et al.|(2018) map. 


For these objects we applied a reddening correction AE(B — V) 
based on the following model 


fi p(r)dr 
fp p@ar’ 


where / and b are the Galactic longitude and latitude, and the 
integral is calculated along the (/, b) direction from the distance 
do at the border of the map and the distance d at which a given 


A E(B — VY, b,d,do) = E(B — V) b) (E.2) 


U http://stilism.obspm.fr 
'8 The map reports precisely this quantity. 


Table E.1. Parameters of the analytical dust's density distribution. Dis- 
tances are in kiloparsecs. 


h. h, 
42 0.088 


T.12Rp 


Yfl Ry Yw 
0.0054 84 0.18 


star is found. The function p is a model of the density distribution 


of the dust. As in (2014), we used 


Ro-R d-z| 
hr kgh; 


p(X) = exp > (E.3) 


where x is the position vector of a star with respect to the Galac- 
tic centre, R the galactocentric distance along the Galactic plane, 
z the orthogonal distance from the Galactic plane, and ky)(R), 
Z,(R) describe, respectively, the flaring and warping of the gas 
disk, 


kg(R) = 1  ygmin(R pi, R — Rei), (E.4) 


Zw(R, $) = yywmin(R,, R- Ry) sin 9. (E.5) 


where ¢ is the Galactocentric azimuth that increases in the direc- 
tion of Galactic rotation with the Sun at $— 0. Table [E.1]reports 
the parameter values. The reddening at infinity E(B — V), in any 
given direction is taken from the Schlegel et al. (1998) map, mul- 


tiplied by the correction factor reported in (2014): 


S(E(B— V)s) = 0.6 + 0.2 03 


1- am (2S === 3 . (E.6) 


The correction described in this section is applied only if the 


reddening calculated up to the|Lallement et al.|(2018) map edge 


was not larger than the reddening E(B-V)., estimated from the 


Schlegel et al.|(1998) map. 


Appendix F: The impact of reddening and 
absorption on the determination of the stellar 
radius 


It is of fundamental importance to estimate the impact of red- 
dening and absorption on the determination of the stellar radius, 
a crucial parameter in the context of transiting planet searches. 
By applying the Stefan-Boltzmann law it follows that the rela- 
tive variation of the apparent stellar radius R, estimated neglect- 
ing the presence of dust, with respect to the true stellar radius Ro 
can be expressed as 


2 
AR (=) 10024+4280) _ 1, (Œ1) 


Ro i Tet 


where Tero is the true effective temperature, and Tep the esti- 
mated temperature of the target star, while A and A BC are the 
absorption and differential bolometric correction in a given pho- 
tometric band, respectively. 

Considering a reddening law, a colour-effective temperature re- 
lation, and an effective temperature-bolometric correction rela- 
tion, it is possible to express the relative radius variation as 
a function of reddening only. For example, adopting the rela- 
tions discussed in Section |5| valid for the Gaia photometric 
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bands, we obtain that, in correspondence to a reddening equal 
to E(Ggp — Grp) = 0.3 mag, the radius of a solar-type star 
would be overestimated by about 3% if both reddening and ab- 
sorption were neglected in the radius calculation, as shown in 
Fig. [ET] sotid line). The individual contributions resulting from 
neglecting either the reddening or the absorption are illustrated 
respectively by the dashed and dotted lines in the same figure. 
Therefore, in correspondence to the same amount of reddening, 
E(Gpp — Grp) = 0.3 mag, neglecting the effect of the reddening 
leads to underestimation of the effective temperature which in 
turn leads to overestimation of the stellar radius by about 3296 
(dashed line). On the other hand, neglecting the effect of absorp- 
tion produces an underestimation of the stellar luminosity, which 
implies an underestimation of the stellar radius by about 2296 
(dotted line). Therefore, reddening and absorption have two op- 
posite effects on the apparent stellar radius which are partially 
compensating each other. 

These results hold for a solar-type star with effective temper- 
ature 7,5—5772 K. It should be kept in mind that this compen- 
sation effect depends on the spectral type of the target source, 
the slope of the adopted colour-effective temperature relation 
(and therefore the adopted photometric system), and the redden- 
ing law. For example, in Fig. [E2] we show that for an M3V star, 
neglecting absorption would produce an underestimation of the 
stellar radius by about 13%, neglecting the reddening would pro- 
duce an overestimation of the stellar radius by about 1196, while 
neglecting both of them would underestimate the radius by about 
496. For a F5V star, neglecting absorption would produce an un- 
derestimation of the stellar radius by about 20%, neglecting the 
reddening would produce an overestimation of the stellar radius 
by about 38%, while neglecting both of them would overestimate 
the radius by about 11%. 

The presence of interstellar matter has a smaller impact on 
M-dwarfs given that these sources are intrinsically redder than 
earlier type stars, and absorption and reddening are more promi- 
nent at shorter wavelengths. Moreover, as we move towards F- 
type stars, the effect of reddening becomes dominant over the ef- 
fect of absorption (because of the strong dependence on effective 
temperature in Eq.[F.1]and the steepening of the colour-effective 
temperature relation, Sect..1) which generally leads to an over- 
estimation of the stellar radius. 
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Fig. F.1. Percentage variation of the apparent radius R estimated ne- 
glecting the presence both of absorption and reddening, with respect 
to the true radius Rọ of a solar-type star (with effective temperature 
Teg—5772 K) as a function of reddening E(Ggp — Gg») (solid line). The 
dashed and dotted lines represent the apparent radius variation resulting 
from neglecting either the reddening (dashed) or the absorption (dot- 
ted), respectively. 
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Fig. F.2. Percentage variation of the apparent radius R estimated ne- 
glecting the presence of both absorption and reddening (dots) for stars 
of different spectral types and for a reddening equal to E(Ggp-Grp)=0.3. 
Squares and triangles represent the percentage radius variation resulting 
from neglecting the reddening or the absorption, respectively. 


